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Abstract 

Maximal lattice free bodies are maximal polytopes without interior 
integral points. Scarf initiated the study of maximal lattice free bodies 
relative to the facet normals in a fixed matrix. In this paper we give 
an efficient algorithm for computing the maximal lattice free bodies 
of an integral (d + l) x d matrix A. An important ingredient is a test 
set for a certain integer program associated with A. This test set may 
be computed using algebraic methods. 

As an application we generalize the Scarf-Shallcross algorithm for 
the three-dimensional Frobenius problem to arbitrary dimension. In 
this context our method is inspired by the novel algorithm by Einstein, 
Lichtblau, Strzebonski and Wagon and the Grobner basis approach by 
Roune. 

1 Introduction 

We will introduce this paper by relating it to the Frobenius problem. Let 
N denote the natural numbers and Z the integers. For G N with 

Zai + ■ ■ ■ + Za„ = Z, the complement of the semigroup Nai + ■ ■ ■ + Na„ in 
N is finite. Its maximum is denoted g{ai, . . . , a„) and called the Frobenius 
number of Consider as an example, the semigroup S generated 

by 6, 10 and 15 in N. The complement of S is 

{1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 14, 17, 19, 23, 27, 29}. 
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and g{6, 10, 15) = 29. Computing the Frobenius number is in general an 
NP-hard problem, which is not known to be in NP (see [11, Appendix A]). 
This phenomenon is perhaps related to the naive method of computing the 
Frobenius number by finding the first consecutive sequence of ai natural 
numbers in the semigroup. The search for a consecutive sequence can be 
simplified through the following classical result due to Brauer and Shockley 
[5]: let Tf he the smallest natural number of the form 

congruent to a given integer / modulo ai, where X2, ■ ■ ■ ,Xn & N. Then 

g{a-i, . . . , a„) = -ai + max{r/ | < / < Oi} (1) 
The reader may find it useful to deduce the classical result 

g{ai, 02) = 0102 - ai - a2 

due to Sylvester (1884) special case. 

In the language of optimization, (1) amounts to solving ai — 1 group 
programs^. These group programs may be interpreted as finding shortest 
paths from in a graph with vertices Z/oiZ and edges suitably weighted 
by 02) • • • jfln- This is basically the graph algorithm in [10] for computing 
the Frobenius number (see also [4]). The graph algorithms have the obvious 
draw-back that they are exponential in the input size in fixed dimension. 

Surprisingly there are polynomial algorithms in fixed dimension related to 
the theory of integral points in convex polytopes. The first such algorithm 
was found by Kannan [8]. Later a remarkable algorithm related to ratio- 
nal generating functions was discovered by Barvinok and Woods [3]. Both 
of these algorithms involve deep mathematical insights, but use quite time 
consuming operations of polynomial complexity in fixed dimension. 

In this paper we present an algorithm for enumerating maximal lattice 
free bodies with facet normals in the rows of a fixed integral matrix. A maxi- 
mal lattice free body is a polytope maximal with respect to having no interior 
integral points. As an application we generalize the Scarf-Shallcross algo- 
rithm [13] for computing the Frobenius number for three natural numbers. 

^The notion of group programs and (graph) algorithms for solving them go back to 
Gomory [7]. 
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Let us briefly recall the beautiful relation between maximal lattice free 
bodies and the Frobenius number: given the coprime natural numbers a — 
(ai, . . . , Qn), Scarf and Shallcross introduced the polytopes 

Kb^{xe R"-^ \Ax<b} 

for b E where A is an n x (n — 1) matrix with columns forming a basis 
of the lattice {v G Z" | av = 0}. The Frobenius number is then given by 

g{ai, . . . , Qn) =max{a6 | b e Z", Kj, maximal lattice free body} 
-J] a,. 

For 3x2 matrices Scarf and Shallcross gave a very efficient algorithm for 
computing the maximal lattice free bodies up to integral translation and 
thereby the Frobenius number for n — 3. 

For integral n x [n — 1) matrices our algorithm generalizes the Scarf- 
Shallcross algorithm. In the context of the Frobenius problem, an algebraic 
version of our algorithm has led to record breaking computations [12]. It is 
possible to extend our algorithm to compute the maximal lattice free bodies 
for arbitrary integral matrices by working with ideal points on facets as in 
[1]. We only treat the simplicial case in this paper. 

This work was originally prompted by an attempt to understand the 
novel algorithm of Einstein, Lichtblau, Strzebonski and Wagon [6] through 
the geometric language in [13]. We thank Daniel Lichtblau and Stan Wagon 
for several discussions related to [6] . 

We also thank Herb Scarf, Bernd Sturmfels and Rekha Thomas for en- 
lightening discussions on neighbors and maximal lattice free bodies. This 
work was carried out at Institute for Mathematics and Its Applications 
(IMA), University of Minnesota. We thank IMA for providing optimal re- 
search conditions during the thematic year on "Apphcations of Algebraic 
Geometry", 2006/7. 

2 Preliminaries 

We begin by recalling a few concepts and a bit of notation from polyhedral 
geometry. First we define the m-dimensional vector max(6i, . . . , 6^) as the 
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coordinate- wise maximum of the vectors bi, . . . ,br G MJ^. For example 




2.1 Polyhedra 

A polyhedron P in is the set of solutions to a finite number of linear 
inequalities in n variables: ajx < hi for G and bi G M, where i = 
l,...,m. We will use the notation PAib) := P = {x E M."- \ Ax < b} 
collecting the normal vectors a\, . . . ,al^ in the rows of an m x n matrix A 
and the bounds bi in a vector b G MJ^. Notice that some of the inequalities 
may be redundant in defining Pa(&)- If 0*2; < bi is redundant we call it an 
inactive facet. 

A polyhedron of the form -Pa(O) is called a polyhedral cone. A polyhedral 
cone is called pointed if it does not contain a line i.e. if Ker(A) = 0. A 
polyhedron, PAib), is called rational if the entries of A and b are rational 
numbers. 

Suppose A is fixed. Then the smallest polyhedron containing vi, . . . ,Vr G 
M"" is given by 

(vi, Vr)A := PA(niax(^f 1, . . . , Avr)). 



3 Integer programs and test sets 

Consider a cost vector c G M". We let IP A,c{b) denote the integer linear 
program 

min{c*x I X G Pa(&) n Z"} = min{c*a; \Ax<b,xe 7/"}, 

where b G and A G M™''". 

Informally a test set for an integer program is a finite set of integral (test) 
vectors, such that a feasible solution is optimal if it cannot be improved 
by moving in the direction of any of the test vectors. Here is the precise 
definition. 

Definition 3.1. A test set for the family IFa,c = {IPa,c(&) I b G W^} of 
integer programs is a finite set T C Z" such that 



4 



(a) c*v < for every v eT. 

(b) fo G PA{b) n Z" is optimal if and only if vq + v ^ PAip) for every v eT. 

Test sets always exist when A is an integral (or rational) matrix. This 
follows for example from the following result due to Cook, Gerards, Schrijver 
and Tardos. 

Theorem 3.2 ([14, §17.3]). Let A he an integral m x n matrix, with all 
suhdeterminants at most A in absolute value, let b be a column m-vector and 
c a row n-vector. Let z be a feasible, but not optimal, solution of maxjca; | 
Ax <b;x integral}. Then there exists a feasible solution z' such that cz' > cz 
and \\ z — z' \\oo< nA. 

For irrational matrices the existence of test sets is more subtle. This issue 
has been addressed in Scarf's theory of neighbors [2]. The following simple 
lemma will be applied in §6.2. 

Lemma 3.3. Let U e GL^(Z). Then T dl/^ is a test set for IP^.c if and 
only if U~^T is a test set for 1^au,cU- 

3.1 Conversion to group programs 

Suppose that A is an integral m x n matrix and b an integral (column) m- 
vector. The integer program IP^ c(^) may be transcribed in the following 
way. A feasible solution a; G Z" to lPa,c(^) corresponds to a vector u G N"^ 
with Ax + u = b. If we let £ denote the subgroup of Z™ generated by 
the columns of A, this means that 1Pa,c(^) niay be formulated as the group 
program Grp^^,{b): 

min {c'Yu 

u = h (mod L) (2) 

u G K\ 

where d G R"* is some cost vector corresponding to c G (recovering x 
from Ax + u — b). In complete analogy with WA,c{b) we have the following 
definition. 

Definition 3.4. A test set for the family {Grp^^ ,,/(&) | b G Z"*} of group 
programs is a finite set T C such that 
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(a) {c'Yv < for every v eT. 

(b) A vector uq G with uq = b (mod £) is optimal if and only if 
for every v e T. 

We use the term group program (as opposed to lattice program) since 
the optimization problem (2) concerns optimizing a function over certain 
representatives of a coset in the group Z"/£. 

In the following we recall how test sets for group programs may be con- 
structed using Hilbert bases for rational cones. This gives in principle an 
algorithm for computing the test set alluded to in Theorem 3.2. 

3.2 Hilbert bases and test sets for group programs 

A pointed polyhedral cone C C M" carries a natural partial order -< given by 

X <y <^=^ y — X & C 

for x,y E C. If C is rational and £ is a subgroup of Z", then the semigroup 
£nC has finitely many -<-minimal elements (see [14, §16.4]). These elements 
are called the Hilbert basis of £ fl C and denoted 3-C(£ fl C). Every element 
of £ n C is a finite non-negative integral linear combination of the elements 
in J{(£nC). 

The Euclidean space decomposes into the 2™ orthants {0^}^™^, which 
are pointed rational cones. For a lattice £ C Z*" we put 

Gr(£) = \J:K{LnOj). 

i=i 

This finite set of vectors is called the Graver basis of £. It is not too difficult 
to prove that {v G Gr(£) | c^v < 0} is a test for the family {Grp^^g(6) |6 e Z'"} 
of group programs. The full Graver basis is a test set for the bigger family 
{Grp£^_p(6) I c e R"*, b e Z"*} of group programs, where both b and c are 
allowed to vary. 



6 



3.3 Generic cost vectors and Grobner bases 

Sufficiently generic cost vectors c e R"* may be viewed as term orders in 
computational algebra. The term sufficiently generic includes the case of 
linearly independent entries over Q in c. In practice we want c^v ^ for 
V E £j inside a sufficiently big ball centered at 0. Usually the cost vector is 
made generic by breaking ties with a term order -<. We will not make this 
precise here but refer the reader to the reference given at the end of this 
subsection. 

The lattice ideal Ic associated to £ is defined as the ideal generated by 
the binomials 

{U""^ -U''' \ vefi} C Q[Ui, ...,Um] 

where v G Z™ is decomposed into vectors v^.v" G with disjoint support 
such that V = — and u"^ denotes the monomial ■ ■ ■ for w = 
{wi, . . . , Wn) G N"*. In this context we have the following result. 

Theorem 3.5. Let c be a sufficiently generic vector in M™. Then a minimal 
Grobner basis of J^, with respect to the term order given by —c is 

-m'^" \veT}, 

where T <Z L is a test set for the family {Gr^i^ J^b) \ b G Z"*} of group 
programs. 

This result is so far the most efficient way of computing test sets for 
integer programs: the algebraic viewpoint makes it possible to apply highly 
optimized algorithms for computing Grobner bases of lattice ideals. A state 
of the art implementation is in the program 4ti2 [http://www.4ti2.de], 
which also contains functions for computing and manipulating test sets for 
integer programs. 

For further information on the relation between Grobner bases and test 
sets we refer the reader to [15]. 

4 Maximal lattice free bodies 

A convex body is a compact convex subset of with non-empty interior. 
We call a convex body lattice free if its interior does not contain any integral 
points. For a beautiful exposition of lattice free bodies we refer to §3 in the 
survey by Lovasz [9]. 
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It is known that a maximal lattice free convex body 5 is a polytope [9, 
Proposition 3.2]. A surprising result due to Bell and Scarf says that B has at 
most 2" facets [14, §16.5]. A very useful characterization of maximal lattice 
free bodies is contained in the following result. 

Proposition 4.1. A polytope is maximal lattice free if and only if its interior 
does not contain any integral points and every facet contains an integral point 
in its relative interior. 

Clearly, Z'^ acts on the set of maximal lattice free bodies in M" by trans- 
lation. We let M{A) denote the set of maximal lattice free bodies of the form 
Pyi(6), where ^4 is an m x n matrix and b e R"* is the right hand side. These 
maximal lattice free bodies are the ones with facet normals in the rows of 
A. For varying b, several of the rows in A may define inactive facets. The 
Z"-action on maximal lattice free bodies restricts to an action on M.{A) by 

z + PA{b)^PA{b + Az). 

It is reasonable to identify maximal lattice free bodies which differ by an 
integral translation. With this convention we have the following result. 

Theorem 4.2. Let A be an integral mxn matrix. Up to integral translation, 
A has finitely many maximal lattice free bodies i. e. 

M(A)/Z" 

is a finite set. 

Proof. Let PA{b) G M(A). We may assume without loss of generality that 
all facets of PA{b) are active. Let ai,...,am denote the rows of A. By 
Proposition 4.1 and integral translation we may also assume that the facet 
corresponding to oi contains in its relative interior. With these reductions, 
we may assume that b — . . . , bmY with 6i = and 6j > for z > 1. 
We now prove that that there can only be finitely many maximal lattice free 
bodies of the form PA{b). 

Let T be the test set in Theorem 3.2 associated with the matrix A. Con- 
sider the integer program 

min{aia; \ Ax <b,x e Z"}, (3) 

Let e > be sufficiently small. Since -Pa(&) has no interior integral points, 
the integer program given by changing the right hand side in (3) according to 
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^2 '■= b2 — e, . . . ,bm '-—bm — ^ has as optimal solution. On the other hand 
is a feasible but not optimal solution for the integer program in (3) when the 
right hand side is changed according to 62 '■= b2 — e, . . . ,bj := bj, . . . ,bm '■ = 
bm — ^-i since Pa(&) contains an integral point in the relative interior of the aj- 
facet. By definition of a test set, there exists Zj e T with aiZj < 0, aiZj < bi 
for j ^ i and ajZj — bj. This shows that a certain subset {^2, ■ ■ ■ , -^r} with 
r — 1 elements of T uniquely defines PA{b) up to integral translation. Since 
there are finitely many such subsets, this proves the claim. □ 



Example 4.3. The 4x2 matrix 



A 



( 1 
-1 


V 



has M(A) = 0. 



\ 


1 



An interesting question is when M(A) ^ 0? This is known to hold if 
A has full rank and every m x m minor is non-vanishing. This condition 
obviously fails in Example 4.3. 



4.1 The simplicial case 

The proof of Theorem 4.2 shows that M(A)/Z" can be computed using a 
test set for the family 

{IPA,a,{b)\beW^} 

of integer programs, where Oi denotes the first row of A. We will describe an 
algorithm for doing this in the more tractable case related to the Frobenius 
problem: assume that A is an integral (rf + 1) x d matrix of full rank such 
that yA = 0, for some y G N^q^. If ao, . . . , are the rows of A and y = 

(l/o, • • • , then yotto H h ydttd = 0. 

These assumptions imply that every d x d minor of A is non-zero and 
every Kj, := PA{b) is either empty, containing one point or a full dimensional 
simplex in W^. The advantage of being in this setting, is that all facets in A 
are active in a maximal lattice free body: let 

Fi{b) ^ {x e Kb\ QiX = bi} 
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denote the facets of a maximal lattice free body Kh, where i = 0, . . . ,d. By 
Proposition 4.1 we may find an integral point in the relative interior of Fo(^)- 
We may therefore assume by translation that is in the relative interior of 
Fo(6). Unless the matrix is sufficiently generic, there may be several choices 
for this translation as the following example shows. 

Example 4.4. The Z^-orbit of the maximal lattice free triangle 



has two representatives containing in the relative interior of -fo(6): 




4.2 Test sets and maximal lattice free bodies 

We are interested in computing a well defined representative of a maximal 
lattice free body in M{A)/Z'^. This means, for example, that we have to 
choose between the two possible candidates in Example 4.4. To make a well 
defined choice of a representative we use the following result. 
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Lemma 4.5. Every maximal lattice free body is an integral translation of a 
unique maximal lattice free body Kf, with the following two properties: 

(a) 6 = (0, 61, ... , bdY e N<^+S where 61, . . . , 6^ > 0. 

(b) is the optimal solution of the integer program 

min{aQ2; | z E U^, a\z < 61 — 1, ... , adZ <bd — 1}, 

where 

Oq := Oq + eai + • • • + e'^a^ 

for e > small, is a perturbation of the first row vector in A. 

Proof. If Kf, is a maximal lattice free body, then the optimal solutions to the 
integer program 

min{ao2; | z E Z'^, 012; < 61 — 1, ... , aaz < &d — 1} 

are the integral points in the relative interior of Fo^b). Therefore they all 
satisfy UqZ ~ 0. By perturbing oq into a'^ for e > small, we identify a 
unique integral solution zq e Fo{b) as the optimal solution to the perturbed 
integer program 

mm{(aQ)z \ z G Z'^, aiz < 61 — 1, ... , adZ < bd — 1}. 

For each maximal lattice free body with in the relative interior, Fq, of -Fo(^) 
this identifies a unique integral point in Fq. The perturbation of ao thereby 
identifies unique representatives of the cosets in M(A)/Z'^. □ 

The key idea is to use the much smaller test set associated with the 
perturbation in constructing representatives of the maximal lattice free 
bodies in M(A): 

Theorem 4.6. Let T be a test set for the family of integer programs given 
by 

min{aQ^ | z G Z*^, aiZ < 61, . . . , ttdZ < bd}, 

where bi, . . . ,bd E Z. If Kb is a maximal lattice free body as in Lemma 4-5, 
then there exists a d-tuple {zi, . . . , Zd) of elements in T such that 

( a) ttiZi > 
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(b) UiZj < ttiZi for j ^ % 

(c) There exists no z & T with a^z < a^z^ for all k — 1, . . . ,d, 

(d) aoZi < 0. 
where i — 1, . . . ,d. 

Proof. This is a straightforward translation of the proof of Theorem 4.2 ex- 
cept for (d): we only have a test set with respect to the perturbed cost vector 
Oq i.e. we have a'^Zi < in (d) but not necessarily aoZi < 0. □ 

4.3 Backtracking maximal lattice free bodies 

In this section we outline a backtracking algorithm for enumerating the 
unique representatives of the maximal lattice free bodies Kj, alluded to in 
Lemma 4.5. Recall that 

Fi{b) = {x e Kh\ GiX = hi} 

denotes the facet given by the i-th row of A. 

By Theorem 4.6, Kf, = (0, Vi, . . . , Vd)A for a certain d-tuple V — {vi, . . . ,Vd) 
of T, such that 

d 

v,eF,{b)\ U FAb). 

Using this observation, there is a simple backtracking algorithm for generat- 
ing these specific tuples discarding most of the d-tuples of T: by definition, 
none of 

{0,Vi)a C ••• C {0,Vi,...,Vd)A 

contain interior points from T. The basic geometric idea is to assign points 
from T to Fi(6), . . . , Fd{b) keeping Proposition 4.1 in mind. Define for v E T 
the set 

Hi{v) = {u E T \ ttiU < ttiv}, 

where 1 < i < d. Assume inductively that we have constructed a (partial) 
tuple V = [vi, . . . , Vi) with the property that 

Vi e Hi{vi) n ■ ■ ■ n Hi_i{vi-i) 

for i > 2 and such that (0,t'i, . . . contains no interior points from T. 

Then we add v E T to V ii and only if 
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(a) V e Hi{vi) n ■ ■ ■ n Hi{vi) 

(b) (0, vi, . . . , Vj, v)a contains no interior points from T. 

If this is not possible, V can never be extended to define a maximal 
lattice free body and we backtrack to V" = {vi, . . . knowing that the 

extension of V" by Vi leads to a dead end. 

By Theorem 4.6, we are sure that B — {0,vi, . . . ,Vd)A is lattice free. 
However if some (aoYvi = 0, B may be contained in a strictly larger lattice 
free body by moving the facet corresponding to Oj. So the backtracking 
algorithm above may generate a superset of the maximal lattice free bodies. 
Note that it is easy to decide from this superset S whether a lattice free body 
is maximal by checking if it is maximal in S. 

5 The reverse lexicographic term order 

Consider the usual lexicographic term order -<iex on vectors in Z" starting 
from the left most coordinate and moving right. Then the reverse lexico- 
graphic term order (for vectors in the subgroup annihilated by a positive 
vector y e N") is given by u ~< v <^=^ —u -<iex —v. 
The perturbation 

u'q = ttQ + eai + ■ ■ ■ + e'^ttd 

considered in Lemma 4.5 has very specific algebraic meaning. The (minimal) 
test set for this perturbation corresponds to a reverse lexicographic Grobner 
basis for the ideal 7^, where £ is the lattice generated by the columns of A. 
This can be read off from the translation (see §3.1) of the integer program 

min{aQ;s | Az < b} 

into the group program 

min{-Mo - eui e'^Ud \ u e N'^+\ m = b (mod £)}. 

When e is infinitesimally small, then 

u <^=^ (-1, -e, -e^, . . . , -e'^){v - u) > 

for u,v e IT-. 
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6 The Scarf-Shallcross algorithm 



In this section we explain how our algorithm specializes to the Scarf-Shallcross 
algorithm for 3 x 2- matrices. First we review the connection between maximal 
lattice free bodies and the Probenius problem given in [13]. 

Consider the Frobcnius problem given by a = (ai, . . . . a„) , where ai, . . . , a^, 
are coprime positive integers. We collect a Z-basis of £ = {f G Z" \ av = 0} 
in the columns of the n x (n — 1) matrix A. The key observation is the 
following. 

Theorem 6.1 (Scarf and Shallcross). Let b G Z". llie integer ah is rep- 
resentable as a non-negative integral linear combination of ai, . . . ,an if and 
only if Kb D Z"-^ = 0, where 

Kb^{xe I Ax < b}. 

Proof. If ab — au for u G N", then b — u G Therefore b — u — Ay ior some 
y G Z"~^ and contains the integral point y. li contains an integral 
point X, then u — b — Ax G N" satisfies au — ab. □ 

Theorem 6.1 shows that 

g{ai, . . . , a„) = max{a6 | 6 G Z" and n Z""^ = 0}. 

In terms of maximal lattice free bodies we have 

51(01, . . . , an) = max{a6 | 6 G Z" and i^^ is a maximal lattice free body } 

6.1 Computing the Frobenius number 

In the notation of the beginning of §6, note that if 

for z G Z"~^ and 6, c G Z", then ac — ab. This imphes that 

giai,.. . ,a„) = max{a6 | Kb + Z""^ G M(A)/Z""^} - 

This means that the algorithm in §4.3 may be used to compute the Frobenius 
number. The output from this algorithm may contain non-maximal lattice 
free bodies along with the representatives of the maximal ones. From the 
point of view of computing the Frobenius number, it is a computational 
advantage to simply compute the maximum above using all these bodies. 
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6.2 n = 3 



In this last subsection we explain how the Scarf- Shallcross algorithm for com- 
puting Probenius numbers in dimension three may be viewed in the frame- 
work of this paper. Wc will argue that the reduction procedure in [13] actu- 
ally computes a test set for a specific integer program. The Scarf-Shallcross 
algorithm is highlighted in Example 6.2 of this paper. To get a feeling for 
the algorithm the reader is encouraged to start with this example. 

Let a — (oi, 02, 03) be three coprime positive integers and put £j — {v & 
\ av = 0}. Then £ has a particularly favorable Z-basis (see [13]) given by 
£. = -|- Zv, where 



u = Aoi and v 





for 7 = gcd(a2,a3) and 7 = Aa2 ~ f^o-s, where A,/x G Z with < /j < 02/7 
and < A < 03/7. By our algorithm we need to compute a test set for the 
family of integer programs given by 

where 61, 62 G Z, to find the maximal lattice free bodies -Pa(&) for the 3x2 
matrix 

/ -7 

A = I Aai —03/7 
\ -(lai 02/7 

The Scarf-Shallcross algorithm ([13], §3) finds a unimodular matrix U (see 
Lemma 3.3) transforming the integer programming problem (4) into 

m.m{c^v I Bv < 6, v G Z^}, 

where 

Ci < 0, C2 < 

Bn > 0, Bu < 
B21 < 0, B22 > 

with B11+B12 > and B21+B22 > 0. For this integer programming problem, 
a test set is {ei, 62, 61+62}. According to Lemma 3.3, {U~^ei, U~^e2, U~^{ei-\- 
62)} is a test set for the original problem. 
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Example 6.2. Suppose we wish to compute g{12, 13, 17) using the Scarf- 
Shallcross algorithm. Then 



-1\ / 

^ = I 12x + 13y + 17z = 0} = Z ( 48 +Z -17 

-36 / \ 13 

The following shows the unimodular transformation U : 




At each step wc add a non-negative of one column to the other preserving 
the sign pattern for the matrix A we wish to reach. We stop when this is not 
possible (as in the last matrix in the example). In this particular example 
we see that every maximal lattice free body is an integral translation of Ki,^ 
or Kh^, where 



bi — max 

and 

62 = max 
The Frobenius number is 







g{12, 13, 17) = max(-12 + 4 • 13 + 1 • 17, -12 + 1 • 13 + 2 • 17) = 57. 
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